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It is shown that the quasiparticle states localized in the vicinity of surface imperfections of atomic 
size can be responsible for the zero-bias tunneling conductance peaks in high-T c superconductors. 
The contribution from these states can be easily separated from other mechanisms using their 
qualitatively different response on an external magnetic field. 
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One of the striking features of high-T c superconductors 
(HTSC) is the presence of zero-bias peaks in the volt- 
age dependence of the tunneling conductance (ZBCP's). 
These peaks, which have been observed in numerous ex- 
periments using low-temperature scanning tunneling mi- 
croscopy (STM) technique ||,f§, planar SIN junctions j§, 
and grain boundaries Josephson junctions |4|, are gener- 
ally considered as a clear signature of the unconventional 
pairing symmetry in cuprates. Thus it is important to 
understand what physical mechanisms can lead to the 
formation of ZBCP's and how the contributions from dif- 
ferent mechanisms can be separated in experiments. 

The most popular model attributes the origin of 
ZBCP's to the Andreev surface bound states (ABS's) ||, 
whose existence in unconventional superconductors is re- 
lated to the fact that quasiparticles reflected from the 
interface see a change in the sign of the order parameter 
along their classical trajectories. The interplay of multi- 
ple Andreev and specular reflections then leads to the for- 
mation of zero-energy bound states in the vicinity of the 
interface. The theory predicts that the ABS contribution 
to ZBCP's is strongly anisotropic reflecting the underly- 
ing symmetry of the order parameter and, in particular, 
is absent for those orientations of the interface for which 
the gap does not change its sign along the quasiparticle's 
trajectory. This result, however, contradicts some of the 
experimental data, in which no significant dependence of 
the ZBCP magnitude on the interface orientation in a- 
b plane has been found Q-q]- On the other hand, the 
analysis of the behaviour of ABS's in an external mag- 
netic field shows that ZBCP's should split symmetrically, 
the splitting being linear in H Q . The experimental sit- 
uation, however, does not lend unambiguous support to 
this prediction. In particular, in Refs. jljj4|, a suppression 
and broadening of the in-plane ZBCP's has been observed 
(see also Ref . Q] , where similar results for the c-axis tun- 
neling were reported). The most plausible explanation 
of the presence of ZBCP's for all surface orientations is 
that the surfaces of real samples are not perfectly flat at 
the atomic scale, so the incident quasiparticles are not 
reflected specularly, but rather get scattered in all di- 
rections resulting in the formation of ABS, albeit with 
a smaller spectral weight [ppl. However, the absence of 



ZBCP splitting in magnetic field cannot be explained in 
the framework of existing theories. 

In this article, we propose a new mechanism for the 
formation of zero-bias anomalies in d-wave supercon- 
ductors which does not rely on the existence of ABS's. 
Briefly, our idea is that a significant contribution to 
ZBCP's, at least for some in-plane surface orientations, 
can come from the states which are localized near atomic- 
scale surface imperfections, and thus are very different 
from ABS's which propagate along the surface. In con- 
trast to the previous approaches, we use an essentially 
non-quasiclassical way of modeling the surface rough- 
ness, namely, we assume that there are strong defects 
at the surface of a lattice superconductor, such as miss- 
ing atoms. It is known that a single scalar impurity has a 
notable effect on the bulk e?-wave superconducting state, 
creating an impurity bound state (IBS) in its vicinity, 
whose energy and width tend to zero in the limit of strong 
impurity potential 0. Experimentally, IBS's manifest 
themselves in the existence of sharp zero-bias peaks in 
the voltage dependence of the differential tunneling con- 
ductance, which have been observed recently in beautiful 
STM experiments on BSCCO compounds |i"0[ . We will 
show that ZBCP's in HTSC can be attributed to the 
formation of zero-energy IBS's by strong surface defects, 
similar to those in the bulk. These states possess the 
two desirable features: (i) they exist for the "anti-node" 
surface orientations, and (ii), in contrast to the Andreev 
states which split in magnetic field, the IBS peaks remain 
centered around the zero energy, but get suppressed and 
broadened. 

Let us consider a single missing atom or a strong re- 
pulsive point-like impurity at a (100) surface of a two- 
dimensional d-wave superconductor. (We consider this 
geometry because it is the simplest case with no ABS's, 
see the discussion in the end of the article.) Because 
of the inherent non-locality of the order parameter in a 
superconductor with higher angular momentum pairing 
and the shortness of the coherence length in HTSC, which 
is typically of the order of lnm, it is convenient to use 
the lattice representation of the Hamiltonian: 



W=£ct(r)tf(r,r')C(r') 
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where C(r) = (cf(r), c|(r)) T are Nambu operators, and 
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is the Bogoliubov-de Gennes (BdG) operator with 
#o(r,r') = -i(r,r') + U(r)5 rr , + fi5 rr ,. The lattice 
sites are labeled by r, and (rr') means nearest neigh- 
bors. We choose the gauge in which the order parame- 
ter is real, so the hopping matrix element has the form 
t(r, r') = £e lp "( r-r The magnetic field is directed along 
the c-axis, and the sample is in the Meissner state, so 
the screening supercurrent decays in the bulk and can 
be assumed to be constant in the vicinity of the surface: 
p s = — (e/c)A(x = 0) = (e/c)HXb, where A is the Lon- 
don penetration depth. The impurity scattering is de- 
scribed by the scalar potential U (r) = uS r o with u — > oo 
(the unitary limit), which describes a surface irregularity 
of atomic size. The mean-field order parameter A(r,r'), 
corresponding to d x 2_ y 2 symmetry, is equal to +Ao for 
r = r' ± a, and — Ao for r = r' ± b, where a, b are unit 
vectors of the square lattice (see Fig. |l|). We do not cal- 
culate the order parameter self-consistently, thus A is 
assumed to be constant. This assumption is valid only 
for the (100) surface orientation, which follows from the 
boundary condition for the order parameter in the ab- 
sence of an impurity . The situation could be further 
complicated if a suppression of the order parameter near 
the impurity site was to be taken into account. Although 
the numerical investigation of the self-consistency equa- 
tions shows that such an effect does exist |l2j , we neglect 
it because it leads only to a renormalization of the effec- 
tive impurity strength towards the unitary limit Jl3| . 

The quantity measured in tunneling experiments is the 
differential conductance, which is proportional to the lo- 
cal density of states (DoS) 



N(r,oj) = — ImGf,(r,r;w) 

7T 



(3) 



where G R is the retarded Gor'kov-Nambu matrix Green's 
function. In the presence of a point-like scalar static im- 
purity, G R = G R + G R TG R , where G R is the Green's 
function of a half-infinite clean superconductor, and the 
T-matrix is given by T(u>) = ut^ [1 — uq(uj)t^\ . Here 
g(ui) = G R (Q,0;ui), and n are the Pauli matrices. The 
poles of the T-matrix correspond to the energies of the 
impurity-induced quasiparticle bound states. A surface 
vacancy (missing atom) prevents electrons from resid- 
ing at its site and thus can be thought of as an in- 
finitely strong repulsive impurity, which corresponds to 
u —t oo. In this limit, the T-matrix has a compact form 
T(lu) — — g _1 (w), so the local DoS induced by a surface 
defect is given by 

SN(r, w) = - Im [G R (r, 0; ^g- 1 (u)G R (0, r; w)] . (4) 
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The unperturbed Green's function G R can be ex- 
pressed in terms of the eigenfunctions ^(r) and eigen- 
values E a of the BdG operator (0): G r (t\, r 2 ; w) = 



J2 a *a(ri)*^(r 2 )/(w + - E a ), where u+ = u + iO. An 
essential ingredient of our theory which makes it different 
from the previous work on IBS's in the bulk, is the neces- 
sity to impose some boundary conditions on the quasipar- 
ticle wave functions at the superconductor- vacuum inter- 
face. If the surface coincides with the x — plane, then 
the boundary conditions are ^f\i(x = —d,y) — 0, where 
d is the lattice constant. It is worth mentioning here that 
if we used a continuum model, then the wave function 
would be required to vanish right at the interface, i.e. at 
x = 0. In this case, a point-like surface impurity would 
not have any effect, so a more complicated approach to 
dealing with surface roughness would be necessary Jl4) . 
In our theory, using the lattice model allows one to incor- 
porate both the non-locality of the order parameter and 
the surface roughness on the atomic scale in a simple and 
natural way. The authors of Ref. |H| used a somewhat 
similar approach in their numerical investigation of the 
self-consistent solution of the BdG equations. Our model 
differs from that of Ref. in several important aspects: 
(i) we clarify the physical reason behind the appearance 
of ZBCP's, namely the formation of bound states, (ii) 
our configuration of surface defects and the way of im- 
posing the boundary conditions at the interface are dif- 
ferent, and (iii) we study the magnetic field response of a 
surface IBS, the problem which has not been addressed 
before. 

The quasiparticle wave functions are considerably 
modified in the presence of a surface, and the Green's 
function satisfying the boundary conditions takes the 
form 

f*/ d dh r n / d dk 

G *(rr,r 2 ;u;)=4/ ^ f*G?(k,«) 

x sink x (x 1 +d)smk x (x 2 +dy k y (vi - V2 \ (5) 



where 



G£(k, W ) 



(w+ - v k p s )r + £ k T3 + A k ri 



(uj+ - v k p s ) 2 - £ 



(6) 



is the Green's function of a bulk clean superconductor 
affected by a "Doppler shift" in the quasiparticle en- 
ergy, £ k = — 2£(cos k x d + cos k y d) — fj, is the normal state 
excitation spectrum, v k — V k £ k is the Fermi velocity, 
and A k = 2Ao(cos/c2;0? — cos k y d) is the momentum- 
dependent superconducting gap. We use the values of the 
parameters typical for YBCO compound: t — 1857716^, 
H = 0.51eV, A = 15meV. 

The energies of the impurity-induced surface bound 
states satisfy the equation det g(iv) = 0, whose solutions 
can be complex. We are interested in the case of strong 
impurity scattering and small supercurrent, so that the 
relevant energies are expected to be small compared to 
the magnitude of the gap. At oj, vpp s <§; A , the momen- 
tum integrals are restricted to small vicinities of the two 



gap nodes (fco, fco) and {ko, — ^o)i where coskod — — ///4t, 
and = 2\/2idsinfco<i is the Fermi velocity at the gap 
nodes. Introducing the notation z = u>/Aq, we have 
g{uj) = (l/4td 2 )F(z)To at complex u, where 

F(z) = —iz H (z + z s ) In (2 + z s ) 

7T 

+-{z - z s )\n(z - z s ) (7) 
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with z s = VfPs/ V^Aq. Two logarithmic branch cuts are 
chosen to go down from z = ±z s parallel to the negative 
imaginary axis. The equation for the spectrum of bound 
states in the unitary limit takes the form F(z) = 0. In the 
absence of supercurrent, F(z) — > Fq(z) = {2/i:)z\nz — iz, 
and the above equation can be easily solved, the solution 
being zq — 0. Both the real and imaginary parts of the 
bare IBS energy vanish, which means that there is a zero- 
energy bound state in the vicinity of a surface vacancy. 
The presence of this state gives rise to a sharp peak in 
the DoS, in full analogy to the situation in the bulk ||. 
At nonzero magnetic field, the dominant energy scale in 
the unitary limit is provided by the Doppler shift VfPs- 
It can be checked that the spectral equation with F(z) 
given by (Q) has only one solution in the complex plane: 
zq = — iirz s /2\ lnz s |, so the IBS energy ujq — zqAq has 
the form 



5N(u, Ps ) = n i ImI(ij)F- 1 ( — 



(9) 
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From Eqs. (||) we see that IBS in the unitary limit is de- 
stroyed by supercurrent and replaced by a resonance peak 
centred around zero energy, whose width depends non- 
analytically on H, proportional to iJ(ln/f) _1 . The phys- 
ical reason for this is clear from Eq. (||): at p s ^ 0, the 
bulk DoS does not vanish at oj — 0, being proportional to 
vfPs, which leads to a much stronger hybridization be- 
tween IBS and the bulk states. It can be shown that at 
large but finite impurity strength, the surface IBS peak 
is shifted away from zero energy, the shift being propor- 
tional to H 2 . However, if VfPs is larger than the bare 
IBS energy Q, the results for the unitary limit are re- 
covered. This behaviour is similar to what is expected 
for impurities in the bulk |l6f| . 

In order to visualize our results and facilitate the com- 
parison with experiments, we compute numerically the 
surface DoS for a finite concentration of surface defects. 
It is easy to see that Nq(0,u>) — > at u — > 00. For this 
reason, in order to study the surface IBS, one should 
calculate either the local DoS, for example at one of 
the nearest neighbors of the impurity site, which can be 
probed by STM technique, or the total interface DoS 
measured in planar tunneling experiments. Here we con- 
centrate on the latter case, which is obtained by putting 
r = (0, y) in Eq. ([!]) followed by the integration over 
y. The resulting contribution to the total DoS from the 
surface IBS's is 



where n, is the linear concentration of surface defects, 
J(w) - (4td 2 /n) fdy [G*{y,w)G*{-y,w)] xl , and F(z) is 
given by Eq. (In). We assume a random distribution 
of defects and neglect the quantum interference effects, 
which gives rise to a prefactor rii on the right-hand side of 
Eq. (g) (see the discussion below). At (lu,vop s )/A — > 0, 
_F _1 (cj) is singular, whereas is not and can be re- 
placed by its value at u — p s = 0, which is real. Wc 
have plotted the results in Fig. ||. The plot confirms 
that there is a sharp zero-bias peak due to IBS's in the 
in-plane tunneling conductance, and the magnetic field 
leads to the broadening of this peak near the zero bias, 
and the suppression of its magnitude. This behaviour is 
in a stark contrast to what is expected for the Andreev 
states. It should be mentioned that the numerical results 
of Ref. [[l7| show that the contribution of bound states 
at randomly distributed strong impurities to ZBCP's in 
the c-axis planar tunneling is either negligible, or leads 
just to a finite conductance at zero bias. We have stud- 
ied a different setup, namely the in-plane tunneling, and 
come to opposite conclusions, a possible explanation be- 
ing that the surface disorder configuration in our system 
favors the appearance of a significant IBS contribution, 
because all the defects lie at the same line - the (100) 
interface. 

We would like to emphasize here that the surface im- 
purity states is not the only mechanism that can lead to 
the observation of ZBCP's for a nominal (100) orienta- 
tion It is the magnetic field response of ZBCP's 
that should help determine which mechanism gives the 
dominant contribution. As said in the introduction, sup- 
pression and/or broadening of the in-plane ZBCP's with- 
out any trace of splitting has been seen in some tunneling 
experiments (see e.g. Ref. Q where the STM results for 
(100) YBCO films were reported), which qualitatively 
agrees with our predictions. However, the possibility of 
a quantitative comparison of the experimental data to 
our results strongly depends on the details of the sam- 
ple preparation and the surface quality. The basic as- 
sumption of our model is that the microscopic surface 
roughness can be described in terms of strong potential 
scatterers of atomic size, which is likely to be the case for 
flat surfaces with missing atoms or steps. On the other 
hand, this assumption is definitely wrong if the charac- 
teristic size of the surface imperfections is larger than 
the coherence length £o- An ideal experimental test of 
our model should be performed on a high-quality (100) 
surface in a broad range of parallel magnetic fields. 

In this article, we concentrated on the zero-energy 
IBS's which are formed in the vicinity of (100) surfaces. 
For other surface orientations, the situation is compli- 
cated by the presence of ABS's, which also have zero en- 
ergy and should therefore experience a strong hybridiza- 
tion with IBS's. The result of this interplay is not clear a 



priori and requires a separate investigation, which is be- 
yond the scope of the present study. Other factors which 
can potentially threaten our results in the presence of a 
finite concentration of defects are the multiple-scattering 
interference effects and the self-consistent order parame- 
ter variation E3|. Although their role certainly deserves 
further analysis, one can always assume that if the sur- 
face disorder concentration is sufficiently small, then the 
characteristic energy scales of our problem, such as the 
Doppler shift, can be made greater than those at which 
the "dangerous" effects mentioned above come into play. 

In conclusion, we have proposed a new mechanism 
of the formation of zero-bias peaks in HTSC. We pre- 
dict that the strong defects at an "anti-node" surface 
can lead to the creation of zero-energy localized states 
in their vicinities, whose properties significantly differ 
from those of the Andreev surface states. These local- 
ized states should manifest themselves by the presence 
of sharp zero-bias peaks in tunneling experiments, which 
get suppressed and broadened in an external magnetic 
field. 
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FIG. 1. Lattice model for a surface defect in a d-wave su- 
perconductor in the presence of magnetic field. 
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FIG. 2. The change in the total interface DoS induced 
by IBS's, as a function of energy at increasing supercurrent 
(rii — 0.1). A sharp peak at zero field is not shown. 



